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ABSTRACT 

We show that the far-IR properties of distant Luminous and Ultraluminous In- 
fraRed Galaxies (LIRGs and ULIRGs) are on average divergent from analogous sources 
in the local Universe. Our analysis is based on Spitzer MIPS and IRAC data of 
L/a>1O 1o L0, 70 /Ltm-selected objects in the 0.1<z<2 redshift range and supported 
by a comparison with the IRAS Bright Galaxy Sample. The majority of the objects 
in our sample are described by Spectral Energy Distributions (SEDs) which peak at 
longer wavelengths than local sources of equivalent total infrared luminosity. This shift 
in SED peak wavelength implies a noticeable change in the dust and/or star- forming 
properties from z~0 to the early Universe, tending towards lower dust temperatures, 
indicative of strong evolution in the cold dust, 'cirrus', component. We show that these 
objects are potentially the missing link between the well-studied local IR-luminous 
galaxies, Spitzer IR populations and SCUBA sources — the z < 1 counterparts of the 
cold z > 1 SubMillimetre Galaxies (SMGs) discovered in blank-field submillimetre sur- 
veys. The Herschel Space Observatory is well placed to fully characterise the nature 
of these objects, as its coverage extends over a major part of the far-IR/submm SED 
for a wide redshift range. 

Key words: galaxies: general galaxies: high-redshift galaxies: starburst galaxies: 
photometry infrared: galaxies 
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1 INTRODUCTION 

In the past 30 years, the field of infrared astronomy has insti- 
gated immense breakthroughs in our understanding of the 
cosmos, at the same time driving pioneering technological 
advances (see review by Low, Rieke & Gehrz 2007). Dis- 
covering and resolving a large part of the Cosmic InfraRed 
Background (CIRB) (Stanev & Franceschini 1998, Gorjian, 
Wright & Chary 1999) led to valuable insights on the na- 
ture of the sources that are responsible for more than half 
the cosmic electromagnetic energy density originating from 
star- formation and Active Galactic Nuclei (AGN). Only 1/3 
of the CIRB is attributed to the local galaxies discovered in 
abundance by the InfraRed Astronomical Satellite (IRAS), 
which implies that a major part of it is composed of red- 
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shifted infrared light originating from galaxies at high red- 
shift. Moreover, the CIRB intensity and peak wavelength 
suggest that the main contributors are dust-obscured pop- 
ulations whose bolometric energy output is dominated by 
emission in the infrared. As studies have repeatedly con- 
firmed, these galaxies hide an immensely active Interstel- 
lar Medium (ISM); they are the ultimate stellar nurseries, 
many also hosting deeply embedded and rapidly growing 
AGN (e.g. Soifer et al. 1986, Sanders et al. 1987, Genzel et 
al. 1998, Tacconi et al. 2002). 

Following these discoveries, a classification scheme 
was introduced based on the total infrared luminosity 
(L/_r) in the 8-1000 /im spectral region and four main 
classes of IR-luminous galaxies were defined: StarBursts 
(SBs, 10 10 <L IR < 10 11 L Q ), Luminous InfraRed Galax- 
ies (LIRGs, 10 11 < Lin < 10 12 L ), UltraLuminous InfraRed 
Galaxies (ULIRGs, 10 12 < L IE < 10 13 L©) and HyperLumi- 
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nous InfraRed Galaxies (HyLIRGs, Lir > 10 13 Lq), e.g. 
Sanders & Mirabel (1996), Genzel et al. (1998). This clas- 
sification is non-arbitrary and a good representation of the 
generic characteristics of these sources, since the transition 
in energetic output from 10 10 to 10 14 Lq is a consequence 
of a change in dust properties (e.g. Noll et al. 2006), star 
formation history (e.g. Cole et al. 2000, Brinchmann & El- 
lis 2000), stellar mass (e.g. Bundy et al 2006), morphology 
(e.g. Murphy et al. 1996, Wang et al. 2006) and the role of 
the central black hole (e.g. Veilleux, Sanders & Kim 1997, 
Brand et al. 2006). 

With data from Spitzer, doors to the high redshift Uni- 
verse were opened, both in terms of resolving the CIRB 
(e.g. Lagache, Puget & Dole 2005, Dole et al. 2006) and 
accurately pinpointing the evolution of the infrared luminos- 
ity function (e.g. Le Floch' et al. 2005, Huynh et al. 2007), 
as well as identifying the key processes responsible for the 
extreme luminosities of these galaxies (e.g. Franceschini et 
al. 2003). In terms of extragalactic science and specifically 
deep Spitzer surveys, the capabilities of MIPS have been 
extensively exploited in the mid-IR, enabling the nature of 
the 24 fj,m population to be characterised up to z~3 (e.g. 
Chary et al. 2004, Marleau et al. 2004, Le Floch et al. 2004; 
2005, Houck et al. 2005, Dole et al. 2006, Marcilac et al. 

2006) . Work on the 160 /im population has been limited to 
source count analysis (Dole et al. 2004a, Frayer et al. 2006a), 
whereas 70 /im-selected sources have also been characterised 
with respect to their IR colours, mid-IR spectra and the in- 
frared luminosity function (Dole et al. 2004a, Frayer et al. 
2006a; 2006b, Huynh et al. 2007, Brand et al. 2008). 

The difference in resolution and sensitivity between the 
three MIPS bands has not enabled a true unification of 24, 
70 and 160 /im populations and hence the range and types 
of galaxies that dominate each one. For example, all Spitzer 
sources detected at 70 /im have 24 jj,m counterparts, but the 
reverse is not true, as the sensitivity of MIPS at 24 fim is 
up to a factor of 100 higher than at 70 fim, considerably 
greater than the f7o/f24 flux density ratio for the average 
IR- luminous galaxy at z>l (typically < 30, Papovich et al. 

2007) . In addition, there is a wide range of sources with weak 
far-IR emission that would form part of a 24 /jm survey, but 
would not necessarily have detections in the MIPS 70 or 
160 /im bands, such as quiescently star-forming spirals, low 
redshift early-type galaxies, optically-selected quasars and 
obscured AGN (e.g. Yan et al. 2004, Desai et al. 2008). On 
the contrary, due to the lower sensitivity of MIPS at 70 /im, 
the majority of far-IR selected objects conform to one main 
type and are thus defined by large amounts of dust, high 
IR to optical ratios and high star formation rates; they are 
predominantly SBs, LIRGs and ULIRGs (Symeonidis et al. 
2007; 2008). This brings a high degree of homogeneity in 
the 70 (im population, clearly benefiting any study of its 
evolution with redshift. 

Infrared Spectral Energy Distributions (SEDs) for 
Spitzer selected sources have been examined with and with- 
out individual 70 and 160 /im detections, the latter via 70 
and 160 ftm stacking on the 24 /im positions (Papovich et 
al. 2007, Zheng et al. 2007, Bavouzet et al. 2008) and the 
former with data from various surveys, such as the Spitzer 
Infrared Nearby Galaxies Survey (SINGS) (Kennicutt et al. 
2003), the Spitzer Wide-area InfraRed Extragalactic sur- 
vey (SWIRE) (Lonsdale 2003; 2004), the All- wavelength Ex- 



tended Groth Strip Survey (AEGIS) (Davis et al. 2007), the 
Bootes field of the NOAO Deep Wide Field Survey (ND- 
WFS) and the Spitzer First Look Survey (FLS) (see Dale 
et al. 2005, Rowan-Robinson et al. 2005, Symeonidis et al. 
2007; 2008, Bavouzet et al. 2008 for work on SEDs). How- 
ever, due to Spitzer's limited sensitivity and angular reso- 
lution at the long wavelengths, the exact shape of the SED 
for various galaxy types has not been examined to the detail 
merited for the high redshift Universe. For non-local sources, 
the short wavelength (Wien) side of the SED peak has been 
mapped in the mid-IR and far-IR up to z~2 and z~l, re- 
spectively, whereas submm studies have examined the long 
wavelength (Rayleigh- Jeans) side up to z~3. Although it 
seems that cumulative data over such a wide wavelength and 
redshift range would have enabled the complete character- 
isation of the IR SED, the different sensitivities of mid-IR, 
far-IR and submm surveys do not favour homogeneneous in- 
vestigations and as a result, many details on the SED are 
missing. Moreover, these 'bolometric' studies do not provide 
detailed information on the key intragalactic environments 
which dominate the energy budget and as a result, funda- 
mental questions with respect to the evolution of the prop- 
erties of obscured galaxies still remain unanswered. 

The natural question that follows is whether local in- 
frared galaxies are the same, similar for the most part or 
instead very different to their high redshift counterparts. 
Here, we approach this topic by comparing the SEDs of lo- 
cal (z~0) infrared galaxies with galaxies at z>0.1, selected 
at 70 /im, based on the preliminary results of Symeonidis et 
al. (2008, hereafter S08). In section 2 we define the sample 
and outline the various selection criteria. Sections 3 and 4 
form the main body of our analysis, where we describe our 
methodology and results, paying particular attention to any 
biases that could have been introduced. In section 5 we dis- 
cuss the implication of our results, making comparisons to 
other studies. Finally, the summary and conclusions are pre- 
sented in section 6. Throughout we use a concordance model 
of Universe expansion, H = 70 kms^Mpc" 1 , Q M = 0.3 
and Q A = 0.7 (Spergel et al. 2003). 



2 SAMPLE SELECTION 

This work is based on Guaranteed Time Observations 
(GTO) using the Multiband Imaging Photometer (MIPS, 
Rieke et al. 2004) on the Spitzer Space Telescope (Werner et 
al. 2004) in the EGS field (~0.5 deg 2 ) (Davis et al. 2007) and 
13 H XMM-Newton/Chandra Deep Field (~0.6deg 2 , centred 
on 13 h 34 m 37 s + 37°54'44", Seymour et al. 2008). We se- 
lect the sample at 70^im and retrieve 178 objects down to 
~4mJy (5a) in the EGS and 244 objects down to ~8mJy 
(5cr) in 13 H Field. As these two surveys cover similar ar- 
eas and share a similar photometric completeness limit at 
f7o ~ lOmJy, we are able to combine the data into one con- 
gruous sample, where 4 < {70 < 290 mjy. 

2.1 Source extraction and photometric coverage 

For details on MIPS data reduction and source extraction, 
we refer the reader to Symeonidis et al. (2007) and Sey- 
mour ct al. (2009), for the EGS and the 13 H field respec- 
tively. Each 70 /um object was matched to the brightest 
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24 /im source within 8", in accordance with the MIPS PSF. 
Given the relative sensitivity limits (~ 110/iJy for 24/im, 
~10mJy for 70/mi), two sources near the 70/tm 5a flux 
density limit with no 24 /jm counterparts, were considered 
spurious and discarded. The 160 /im photometry was ex- 
tracted with the IRAF /DIGIPHOT package on the known 
positions of the 70 /im sources, enabling us to go nearer the 
confusion limit, estimated to be on the order of ~40mJy 
at 160 /im (e.g. Dole, Lagache & Puget 2003). We used an 
aperture radius of 40 "(equivalent to the 160 /im diffraction 
limit) and the aperture corrections specified by the MIPS 
team (1.884 for a 40"radius and 40-75"background annulus; 
|http:/ /ssc. spitzer.ca ltech.edu/mips/apercorr/ ). We tested 
our method of source extraction on the Spitzer Cosmic Evo- 
lution Survey (S-COSMOS, Frayer et al. 2009) 160/tm im- 
age, with the known positions of the S-COSMOS 160/im 
sources and found our flux density estimates to be consis- 
tent, within 10 per cent (r.m.s.) of the S-COSMOS MOPEX 
(Makovoz & Marleau 2005) estimates. 

The MIPS photometry was supplemented with In- 
fraRed Array Camera (IRAC) 8 fim data to ensure a more 
extensive coverage of the infrared SED. The overlapping area 
of the 8, 24 and 160 /jm images adds up to about 65 per cent 
of the ECS 70 /im survey and 54 per cent of the 13 H field 
70 /im survey. Requiring IRAC coverage results in a subset 
of 114 70 fim sources in the ECS and 132 in the 13 H field, 
instead of the 178+224 described earlier. 

2.2 Redshifts for the 13 H field 

Photometric redshifts were estimated for the 13 H field, using 
multi-band deep optical and near-IR imaging: u* , B, g', R, 
i', I, z' , Z, J, H, and K, from various ground based obser- 
vatories (CFHT, Subaru, INT and UKIRT) and 3.6, 4.5, 5.8, 
8.0/im from Spitzer's IRAC. The requirement of area cov- 
erage in at least 4 optical/near-IR bands, reduced the num- 
ber of available 70 /im sources from 132 to 108. The cross- 
matching was done using the 24 /im positions each 70 /tm 
source was identified with. Photometric redshifts were es- 
timated with the HyperZ (Bolzonella et al. 2000) template 
fitting code, including the set of galaxy and AGN SED tem- 
plates from Rowan- Robinson et al. (2008) and setting optical 
extinction as a free variable. All aperture magnitudes were 
also corrected for aperture losses. Moderate optical extinc- 
tion was permitted for the AGN and late type galaxy tem- 
plates, and heavy extinction (up to Ay = 10) was allowed 
for the starburst templates. For a full description of the pho- 
tomety extraction and redshift fitting process see Dwelly et 
al. (in prep.). 

Within the final overlap area of the Spitzer MIPS /IRAC 
and ground-based optical/near-IR surveys, 39 good quality 
spectroscopic redshifts are also available, as part of an on- 
going program to identify the optical counterparts of faint 
X-ray and radio sources in the 13 H field, using various long- 
slit and multi-object spectrographs at the WHT, CFHT and 
Keck telescopes. The reliability of the photometric redshifts 
was checked against the spectroscopic redshifts, finding high 
consistency with an r.m.s. scatter in redshift of only ~ 0.05. 

Finally, as the aim of this work is to evaluate the proper- 
ties of non-local 70 /im selected sources, we made a redshift 
cut at z=0.1, rejecting the 26 low redshift resolved bright 
galaxies in the 13 H field. The final 13 H field redshift sample 
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Figure 1. The combined redshift distribution of the 70/tm sample 
(zmeun ~ 0.5, Zmadian ~ 0.42; solid black line) from the individual 
distributions of the EGS population (dot-dashed red line) and the 
13 field population (dashed blue line). The difference between 
the two redshift distributions is likely due to sample variance. 

consists of 82 sources, 34 with spectroscopic redshifts and 
48 with photometric redshifts. 

2.3 Redshifts for the EGS 

Spectroscopic redshifts for the EGS were acquired by the 
Deep Extragalactic Evolutionary Probe 2 (DEEP 2) survey 
(Davis et al. 2003). DEEP 2 targeted sources in the EGS 
with a magnitude limit of Ras = 24.1. The brightest, most 
extended sources were weeded out of the spectroscopic tar- 
geting algorithms, as their photometry was suspect and they 
could potentially saturate the detector. Thereafter, the re- 
trieval of redshifts was random, at ~60% sampling and at 
the time of writing there are 66 good quality redshifts avail- 
able for the 114 sources in the MIPS/IRAC overlap areas. 
These are in the range 0.1 <z< 1.9, with a mean of 0.55 
and a median of 0.48. The exclusion of the brightest R<18 
sources implies that there is an under-representation of ex- 
tended, low-redshift galaxies. However, as we apply a red- 
shift cut at z=0.1 (see section |2~I2")) . these would have been 
eliminated from our final sample anyway. On the faint end, 
there are only 4 objects fainter than R=24.1 in the full sam- 
ple, i.e. only two would have made it in the redshift sample, 
so overall there is high consistency between the two samples. 

The limited optical photometry for the EGS 70 /tm ob- 
jects did not allow us to calculate reliable photometric red- 
shifts and as a result, we only use these 66 sources with spec- 
troscopic redshifts in subsequent work. However, as outlined 
above, this redshift sample is a good represenation of the the 
z>0.1 parent EGS sample. 

2.4 The final 70 /im sample 

All 148 sources, assembled after the various area and red- 
shift cuts described in earlier sections, are detected at 8 and 
24 /tm. Eight sources have signal to noise < 1 at 160 /tm and 
are therefore excluded from subsequent analysis. From the 
remaining 160 /tm measurements, ~80 per cent are >3<r de- 
tections, ~ 62 per cent >5<r and ~47 per cent >7a. 
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The final sample consists of 140 galaxies in the 0.1- 
2 redshift range (z mean ~ 0.5, Zmedian ~ 0.42, hereafter the 
70\im sample). Figure [1] shows the individual and combined 
redshift distributions, the differences likely arising due to 
sample variance. 

Completeness is not a strict requirement for our anal- 
ysis. However the selection criteria we apply, ensure that 
the 70 [im sample is representative of the parent 70 /im pop- 
ulation. Nonetheless, given the strong correlation between 
R-band magnitude and redshift, we examine the possibility 
that the requirement of (optically derived) redshifts modi- 
fies the type of 70 /j,m sources present in the final sample. 
We find that the redshift at which each source becomes too 
faint for inclusion in the 70 /im sample is much lower than 
the redshift at which it drops out of the R-band survey. This 
is to be expected as the R-band limit (< 1/iJy), is 4-5 or- 
ders of magnitude lower than the 70 /im flux density limit 
(~ lOmJy), but at the same time ~95 per cent of the ob- 
jects have fVo/fii <5000. Under these circumstances, sources 
would drop out of the 70 fim selection before dropping out 
of the optical selection, confirming that the sample of 140 
sources we focus on, is representative of z>0.1 70 /im de- 
tected galaxies. 



2.5 The local sample 

In order to compare local and high redshift sources in a 
way that is as unbiased as possible, we evaluate our results 
against 259 galaxies from the IRAS revised Bright Galaxy 
Sample (BGS) selected at 60 /im (hereafter, the local sam- 
ple) by using available photometry from Sanders et al. (2003, 
hereafter S03). The IRAS BGS is particularly suitable for 
such a comparison because it includes all nearby LIRGS 
and ULIRGs and the 60 /im selection criterion is compara- 
ble to 70 /im at z~0.1 where our sample begins. We do not 
weed out any sources which potentially have a significant 
AGN contribution from either the local or the 70 (im sample. 
Apart from the fact that there is a low AGN incidence in far- 
IR selected samples, both populations are defined by pow- 
erful far-IR emission, hence the presence of large amounts 
of dust, explicitly setting star-formation to be a very a sig- 
nificant part of the galaxies' energetics, if not the domi- 
nant. The IRAS BGS consists of 629 galaxies (z<0.08), 
the majority of which are starbursts, but with a significant 
LIRG population of ~ 25 per cent and a very small fraction 
of ULIRGs (~ 2 per cent). We only use galaxies which are 
not flagged as having large uncertainties in the photometry: 
all the available LIRGs and ULIRGs (143 and 18) and 98 
randomly-selected starbursts. As our main conclusions rest 
on the comparison between local and high redshift LIRGs 
and ULIRGs and much less so on SBs, we believe it unnec- 
essary to include all of the IRAS SBs in the local sample. 
Consequently, we rely on a randomly-selected subset, which 
is about 3-4 times the number of starbursts in the 70\im 
sample, to represent the BGS SBs. 

3 SED FITTING 

The SED shape of an IR-luminous galaxy in the 5-1000 ^m 
range (rest wavelength), is dominated by emission from dust, 



with a range of temperatures spanning ~ 3 orders of mag- 
nitude from ~ 1500 K at the shortest wavelengths to ~ 10 K 
at the longest. In the near-IR (~2-5/im), light from stel- 
lar photospheres dominates the continuum, although there 
can be a small contribution from dust (e.g. Lu et al. 2003), 
subsequently becoming negligible shortward of 2 /xm, as the 
maximum grain sublimation temperature is of the order of 
~ 1500 K. The mid to far infrared part of the SED has two 
major contributions: emission from stochastically heated 
Very Small Grains (VSGs), with T>60K, responsible for 
the continuum up to about A ~ 40 ^m and emission from 
large grains in thermal equilibrium, with T< 50 K, respon- 
sible for emission longward of A ~40 ^m (e.g. Efstathiou & 
Rowan-Robinson 1990; 1995, Granato & Danese 1994, Silva 
et al. 1998, Klaas et al. 2001). Accordingly, the rest-frame 
SED peak, which defines the wavelength of maximum en- 
ergetic output (in vi u ), is in the 40-140 /xm range and the 
SED is often approximated by a modified black body (B^v 13 ) 
of T~ 20-60 K and emissivity /3 ~ 1-2. As the bulk of the 
dust mass is in equilibrium, such an approximation does 
represent the average SED shape, however a colder, 'cirrus' 
component associated with the properties and distribution 
of dust in the ISM and the intensity of the diffuse radiation 
field is often missed. If the cold component is significant, yet 
ignored, the integrated energy output of the galaxy will be 
underestimated, but more importantly it will lead to erro- 
neous estimates for physical quantities, such as dust mass, 
emissivity, opacity and temperature — e.g. see Silva et al. 
(1998), Granato et al. (2000) and Siebenmorgen, Krugel & 
Laurejis (2001) where the importance of both the starburst 
and cirrus components is emphasized. 

In this paper, we broadly categorise the sources as 'cold' 
or 'warm', depending on the average temperature that de- 
fines the bulk of the dust mass, which, as mentioned ear- 
lier, is mostly found to be in thermal equilibrium. As a 
result the average dust temperature scales with the in- 
verse of the SED peak wavelength; a cold galaxy SED is 
expected to peak at longer wavelengths, whereas a warm 
SED will peak at shorter wavelengths. Note that this is 
not the same as the classical IRAS 'warm/cold' terminol- 
ogy parametrised as the mid-IR /2s //eo continuum slope. 
The IRAS 'warm/cold' /25//60 criterion is based on the ex- 
pectation that dust heated by an AGN would reach higher 
temperatures and hence lower the continuum slope between 
25 and 60 /mi (or 24 and 70 /mi in Spitzer's case) (e.g. Far- 
rah et al. 2005; Verma et al. 2005, Frayer et al. 2006); it 
has been applied extensively in order to quantify relative 
AGN/starburst contributions in ULIRGs (e.g. de Grijp et al. 
1985, Miley, Neugebauer & Soifer 1985, Soifer et al. 1989). 
What we refer to as cold/warm, is not related to the mid- 
IR continuum slope, but solely to the far-IR properties and 
more specifically, the position of the SED peak. 

3.1 Methodology 

We fit the available IR data for the 70fim and local sam- 
ple with the Siebenmorgen & Krugel (2007, hereafter SK07) 
theoretical SED templates. In S08, we found that these mod- 
els are optimum for fitting a wide range of source types, 
as their formulation employs a configuration of stars and 
dust distributed throughout the volume. In the SK07 for- 
mulation, mid-IR emission originating in dust regions en- 
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Figure 2. The 60/450 vs 100/850 flux density ratio for local 
galaxies. Open diamonds: predicted colours for the local sam- 
ple using the SK07 fits (black for SBs, blue for LIRGs and red 
for ULIRGs). Filled black circles: observations from Dunne & 
Eales (2001, DE01). Note the good agreement between predicted 
and observed colours, verifying both the accuracy of our fitting 
method and the validity of the SK07 models. 



veloping OB stars is deconvolved from the cirrus emission 
due to the general stellar radiation field, allowing the far- 
IR part of the SED to be modeled as a two-temperature 
component; see also Klaas et al. (2001), Dunne and Eales 
(2001) for the benefits of a 2-temperature far-IR formula- 
tion. The SK07 library consists of ~ 7000 templates defined 
by the physically acceptable combinations of the following 
5 parameters: radius of dust emitting region (0.35, 1 and 3 
kpc), visual extinction (Ay = 2-150), total infrared luminos- 
ity (10 10 -10 14 7 L Q ), hot spot dust density (100-10000 cm" 3 ) 
and L 3tars ,o.B/L S i Qrs , tot (40, 60 and 90 per cent). Hot spot 
dust density refers to the density of the dust region envelop- 
ing OB stars. L star3 ,os/L st ar 3 ,tot, represents the percent- 
age of the total infrared luminosity that originates from OB 
stars, compared to the contribution from the general stellar 
population. 

We perform a \ 2 fit on the photometry (8-160 /im for 
the 70/j.m sample; 12-100 fim for the local sample), using the 
entire range of SK07 templates, folding in the photometric 
redshift uncertainties (where applicable). For each object, 
we determine one interesting parameter, namely the wave- 
length of the SED (yLu) peak, calculated using the mean 
(and la uncertainty) of all templates with X; ^ Xmin + 1| 



where i is the template index and Xmin is the minimum \ 2 
value. In the same way we also get estimates for the total in- 
frared luminosity (Lir, 8-1000 /im) and submm flux density 
at 350 and 850 fjm. 

We find that the 70 fim sample consists of 30 SBs, 79 
LIRGs, 29 ULIRGs and 3 HyLIRGs, indicating that down 
to the ~10mJy level, L>1O 1O L0 sources are the sole con- 
tributors to z>0.1 70 /im populations. 

Note that any comparison between IR galaxies at var- 
ious redshifts will only be fair if each luminosity class is 
examined separately. As outlined in section [T] the transition 
in luminosity closely follows changes in the galaxies' phys- 
ical properties, star- formation history, etc., which must be 
reflected in all comparative evaluations, something that we 
maintain throughout this work. 




redshift 

Figure 3. Top panel: The distribution in 60 /im flux density 
for the local sample (black histogram), compared with the dis- 
tribution that the 70 (im sample would have if it were at the 
average redshift of the local sample (z~0.012; ~50Mpc) (green 
histogram). Note that apart from a few low luminosity SBs 
(L<5xl0 10 L ), the 70 um sample would be easily detectable 
within the 5.24 Jy IRAS flux density limit. Lower panel: Total 
infrared luminosity plotted as a function of the maximum redshift 
at which the types of sources that form part of the local sample 
would be detectable within the 70 fim 10 mjy flux density limit 
- SBs in black, LIRGs in blue and ULIRGs in red. Comparing 
with the average Lir and redshift for each luminosity class of 
the 70(im sample (filled circles), reveals that the 70 /im selection 
is quite sensitive to the LIRGs and ULIRGs of the local sample; 
at the high luminosity end, the sources in the local sample can 
be detected to higher redshifts than the average redshift of the 
70 fim galaxies. 



3.2 Selection Issues 

3.2.1 The SED models and fitting accuracy 

We chose the SK07 models, as their theoretical formulation 
allows the representation of a wide range of systems, min- 
imising assumptions on the types of sources that exist. This 
is in contrast to empirical sets of templates, e.g. Chary & El- 
baz (2001), Dale & Helou (2002), which are tailored to local 
galaxies and could potentially predispose the results towards 
the types of SEDs found locally. In S08 we show that this 
is indeed the case and find them un-representative of the 
SEDs of 70 /xm-detected galaxies, especially in the far-IR. In 
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Figure 4. The 160/70 ^m flux density ratio versus redshift for the SBs, LIRGs and ULIRGs of the 70fim sample (filled cirles). The 
dotted line tracks are the Chary & Elbaz (2001) templates and the solid line tracks are the Rieke et al. (2008) templates with luminosities 
of log Ljr = 10, 11, 12. All templates have been convolved with the MIPS 70 and 160 fim filter response curves, since MIPS photometry 
is calculated from the integrated flux in each band and is therefore not strictly monochromatic. The colour coding for both data points 
and tracks is the same: black for starbursts, blue for LIRGs, red for ULIRGs. The errors have been estimated by combining calibration, 
photometric and statistical uncertainties. Note the large scatter in fieo/f70 spanning almost an order of magnitude within each luminosity 
class. Moreover, the fi6o /iyo values for many 70 ^tm-detected LIRGs and ULIRGs are offset by at least a factor of 2 and up to 10 in 
extreme cases from local tracks of the same luminosity, suggesting that these sources are on average colder than their local counterparts. 



contrast, we find that the SK07 models are representative 
of the SEDs of both the local and 70 fim samples. 

To quantify this further, we compare to submm data 
from Dunne et al. (2000, hereafter D00) and Dunne & Eales 
(2001, hereafter DE01), who describe the submm properties 
of local galaxies in the SCUBA Local Universe Galaxy Sur- 
vey (SLUGS) and report 850 and 450 fim observations for a 
fraction of the IRAS BGS. Using the SK07 models we pre- 
dict 60, 100, 450 and 850 /im flux densities for the sources 
in the local sample and evaluate the results against the ob- 
served DE01 far-IR/submm colours (60/450 vs 100/850, fig- 
ure[2]). We find good agreement between our predictions and 
the DE01 observations, which confirms both the accuracy of 
our fitting method and the validity of the SK07 templates. 

3.2.2 The detection limits of the 60 fim IRAS survey and 
70 \im Spitzer survey 

Selecting the sample at a wavelength shortward of the peak 
(the Wien side of the SED), by definition, does not pre- 
dispose the selection towards cold sources, in the same way 



that selecting longward of the peak (the Rayleigh- Jeans side) 
does not predispose towards warm sources. Our 70 fim se- 
lection is well into the Wien side of the SED (for z > 0.2 it 
corresponds to X re st < 58/im), so any bias would favour the 
selection of warmer SEDs, especially for galaxies at z > 0.5. 
The 60 fim IRAS selection probes longer wavelengths in lo- 
cal galaxies than the 70 /jm MIPS photometry for sources at 
z>0.2, so, if anything, it should be more sensitive to rela- 
tively cold SEDs. Also, note that the 160 /jm photometry co- 
incides with the IRAS 100 fim band for sources at 0.4< z <1; 
this applies to all ULIRGs and the higher luminosity LIRGs 
in the 70 fim sample. 

We estimate the 60 fim flux density the 70 fim sources 
would have if at the average redshift of the local sample 
(z ~ 0.012, ~50 Mpc) and compare with the f6o distribution 
for the local sample, finding considerable overlap (figure [3j 
top panel) . It is only the lowest luminosity 70 /xm-detected 
SBs (L< 5x 10 10 Lq) that would not be recovered within the 
IRAS 60 fim flux density limit of 5.24 Jy. This is not surpis- 
ing, as at 70 fim we detect them at the mjy level and at low 
redshifts (z<0.2). Using a (conservative) 70 fim lOmJy flux 
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Figure 5. The relationship between SED peak and total infrared luminosity for the local sample (black filled circles) and 70fim sample 
(red open squares), split into 3 luminosity bins (SBs, LIRGs and ULIRGs), outlined by the vertical dashed lines. The solid line is the 
peak-Ljfl relationship of the Rieke et al. (2008) templates. Note that for a given Lir, many high redshift galaxies have on average 
significantly colder SEDs that peak at longer wavelengths. 



density limit and the SK07 templates matched to the IRAS 
galaxies, we calculate the maximum redshift at which each 
source in the local sample would still be bright enough to 
appear in the 70 /im survey (figure |3j lower panel). As ex- 
pected, the types of sources that make up the local sample 
are easily recoverable at f7o > lOmJy. In fact, at the high lu- 
minosity end (log L_f_R>11.5), the sources in the local sample 
are detectable to higher redshifts than the average redshift 
of the galaxies in the 70/j.m sample. The fact that the en- 
tire range of template matches for L>5xlO lo L0 sources in 
both samples is detectable with the flux density limits of the 
Spitzer/MIPS and IRAS surveys, confirms that the samples 
are well-matched and that our results and conclusions (sec- 
tions 2] and [SJ are not affected by the initial selection. 

The final point to address is the dependency of the IR 
luminosity and position of SED peak on the absolute value 
of the 160 fim flux density. Our detection levels are at least 
5 a at 70 /im and typically 4-5 a at 160 /im, so to simulate 
combined uncertainties we evaluate the effects of a 30 per 
cent change in fi6o for all objects in the 70/im sample. We 
find that the change in luminosity does not exceed 0.2 dex 
and the average peak wavelength shift is less than 10 /im, 
well within the typical uncertainty (see section 14. 2[) . This 



confirms that the absolute value of fi6o does not drive our 
estimates for X pea k and Lj_r. 



4 RESULTS 

4.1 The 160/70/im flux density ratio 

As a first approach we examine the observed f i6o / fVo colours 
of the 70fim sample as a function of redshift (figure [4}, re- 
lying on SED templates to represent the IR colours of local 
galaxies. Accordingly, we chose the Chary & Elbaz (2001, 
hereafter CE01) and Rieke et al. (2008, hereafter R08) em- 
pirical templates, well-representative of the average proper- 
ties of local galaxies, but with some key differences in the 
continuum slope shortward of the peak; R08 consider their 
templates more accurate for high luminosity IR galaxies, 
such as Arp220, as they account for the steeper continuum 
slope seen in the spectra from the InfraRed Spectrograph 
(IRS) on Spitzer. 

The fi6o/f70 ratio directly scales with the position of 
the peak, a relationship which is much clearer once K- 
corrections have been applied (see S08). Here, in order 
to minimise assumptions, we only compare the observed 
fi6o/f70 colour against the locally-derived tracks. However, 
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Figure 6. The distribution in SED peak for the 3 main lumi- 
nosity classes (SBs, LIRGs and ULIRGs). The dotted red line 
represents the local sample and the solid black line represents the 
70/im sample. Although there is overlap in the distributions, a 
significant fraction of high redshift LIRGs and ULIRGs peak at 
longer wavelengths than the local galaxies of the same luminosity. 



by splitting the 70/j.m sample into luminosity classes we can 
evaluate fi6o/f70 more consistently and in a relatively nar- 
row redshift range. The scatter in SB fi6o/f70 is in broad 
agreement with the local SB SEDs, with a similar number 
of sources lying above and below the region defined by the 
tracks. On the contrary, the fi6o/f70 for many LIRGs and 
ULIRGs is displaced from the local LIRG and ULIRG tracks 
by at least a factor of 2 and up to 10 in extreme cases. This 
suggests that a large fraction of these galaxies are colder 
than their local equivalents, as the energetic output at long 
wavelengths is increased relative to the short wavelengths. 
Overall, the fi6o/fVo colour does not seem to change with 
redshift and, although this would be mis-interpreted as a 
no-evolution effect, it is in fact quite the opposite. Examin- 



ing the variation of f i6o / fVo with redshift will invariably have 
a luminosity and if-correction bias folded into it, as at each 
redshift range the survey probes galaxies at different rest- 
frame wavelength and of different L/_r. In fact strong evolu- 
tion is evident, as the high redshift LIRGs and ULIRGs have 
fi6o/f70 ratios comparable to those of local low luminosity 
starbursts as opposed to local LIRGs and ULIRGs. 



4.2 The SED peak 

In figure [5] the variation of SED peak (rest-frame wave- 
length) is plotted against total infrared luminosity for the 
local and 70 /im sample. Broadly speaking, the peak repre- 
sents an average temperature for the bulk of the dust. In the 
nearby Universe, more luminous galaxies are characterised 
by emission from dust at higher temperatures which shifts 
the SED peak to shorter wavelengths, explicitly shown in 
the SED peak-L/_R relationship of the local sample; see also 
Sanders & Mirabel (1996) and the peak-L/_R variation in the 
R08 SED templates, overplotted for comparison. 

The local and high redshift sources, especially objects 
that fall in the LIRG and ULIRG luminosity class, occupy 
almost distinct regions of this plot, signifying that for a given 
L/_R, high redshift galaxies have on average colder SEDs that 
peak at longer wavelengths. Since the SED peak seems to 
evolve from the local to the high redshift Universe, we might 
also expect it to evolve with redshift within the 70\im sam- 
ple, i.e. from z^O.l to ~2. However, we are unable to ex- 
amine this because, as mentioned earlier, the L/_r-z plane is 
not sufficiently sampled. Under the circumstances, we rely 
on the more unbiased comparison between local and high 
redshift sources by examining the distribution in peak wave- 
length separately for the 3 luminosity classes — figure [6] 
SBs in both samples seem to be broadly consistent, as the 
SB population in the 70\xm sample is at quite a low mean 
redshift (z~0.2) and cosmic evolution in the z~0 to z~0.2 
interval is not expected to be prominent. However, there 
are a few exceptions, with some 70 pm-selected SBs being 
colder than local equivalents, but since these sources are 
at low redshifts, this could be a selection effect due to the 
much deeper MIPS survey. As described in section 13.2.21 
this argument does not apply to the higher redshift, higher 
luminosity sources and hence any differences in the types of 
objects picked up by the two surveys will be solely due to 
evolution. 

Figures [5] and [6] show that, although there is some over- 
lap in the types of objects selected in the two surveys, it 
seems that the SEDs of a substantial number of high red- 
shift LIRGs and ULIRGs peak at longer wavelengths than 
their local counterparts. Such a hypothesis was proposed in 
S08, where fitting locally-derived SED templates underesti- 
mated the far-IR emission of a large fraction of high redshift 
objects. Also, it seems that the galaxies in the 70\im sample 
span approximately the same range in SED peak, whatever 
their luminosity class, in contrast to the narrower local dis- 
tribution, which notably shifts towards shorter wavelengths 
at higher luminosities. Although the majority of local LIRGs 
peak at ~75/im, with a smaller fraction at ~100 /im (mean 
^ P eafc~77/im), high redshift LIRGs span a wide range in 
SED peak up to 160 /im, with a mean of ~105^m. This 
effect is even more pronounced in the ULIRGs, where the 
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average peak in the SEDs shifts from an average local value 
of ~60/im to a high-redshift value of ~92 /im. 



4.3 Submm colours 

In sections l4.ll and r4.2l we presented strong evidence that the 
70jim sample is not simply a redshifted version of sources 
in the local Universe. What naturally follows is a compari- 
son with the SubMillimetre Galaxy (SMG) population, com- 
prised mostly of cold galaxies at high redshift. Due to the 
negative 7f-correction and sensitivity limits, SMGs discov- 
ered in blank field surveys are mostly in the 1.5<z<3.5 red- 
shift range (e.g. Chapman et al. 2005), as the 350-850/im 
wavelength range corresponds to 140-340^im < \ res t < 78- 
188/xrn, effectively sampling near the SED peak for a wide 
range of sources. 

Using the SK07 fits, we predict 350 and 850/im flux den- 
sities for the LIRGs and ULIRGs in the 70[im sample, with 
uncertainties derived from the range of best-matched tem- 
plates (see sect ion l3~Tjl and plot the evolution of the faso/fsso 
colour with redshift (figure [7]). We compare with LIRG and 
ULIRG SMG data from Kovacs et al. (2006, hereafter K06) 
and Coppin et al. 2008 (hereafter, C08). Also plotted are 
the R08 tracks, to represent the local IR galaxies. As the 
LIRGs in the 70[im sample are at an average redshift of 
0.4, we predict the sub-mm flux densities they would have 
if their redshifts were increased by 0.7 and 1.5, in order to 
bring them to a similar redshift range as the K06 and C08 
LIRGs (top panel) . Similarly, we calculate what submm flux 
densities we would expect for the ULIRGs (average z~0.8) 
if their redshift were increased by 1.3 and 2 (lower panel). 

The majority of K06 and C08 LIRGs and ULIRGs 
fall below the corresponding R08 tracks, implying a lower 
f35o/f850 ratio; for a given L IR , sources with lower faso/fsso 
are colder, as there is a shift in far-IR emission, with more 
flux coming out at longer wavelengths. A large number of 
sources from the 70 \xm sample have similarly low faso/fsso, 
again mostly below what is denoted by the R08 tracks, in- 
dicative of more cold dust emission, compared to what is 
seen locally. The agreement in fsso/fsso between the SMGs 
and MIPS sources is a key result; it strongly suggests that 
the galaxies we detect at 70/im are the z<l equivalents of 
the cold z>l sources detected in submm surveys. 



5 DISCUSSION 

5.1 Cold Galaxies at z < 1 

We find evidence for evolution in the dust/star-formation 
properties with redshift, initially apparent through the 
fi60 / fro colour and further quantified by the rest frame wave- 
length of the SED peak. More specifically, we show that a 
population of LIRGs and ULIRGs at 0.1<z<2 appears to 
have a cold-dust-associated far-IR excess not seen locally. 
Although we do identify many sources consistent with the 
kind of objects detected in the local Universe, such as the 
well-studied M82 and Arp220 types, our results establish 
a shift in the rest-frame far-infrared SEDs of high redshift 
sources towards longer wavelengths. A simple description for 
this, is that the shapes of the far-infrared SEDs of many high 



redshift sources may resemble more closely those of lower- 
luminosity galaxies locally. 

The strong cold component identified in many of our 
70 /im objects must relate closely to conditions in the ISM, 
as it represents an increase in far-IR/submm flux com- 
pared to local galaxies and hence must represent a change 
in dust properties and/or star formation efficiency. For a 
given total infrared luminosity, sources that are described 
by colder SEDs may consequently be characterised by larger 
dust masses and/or higher dust opacity, an increase in grain 
emissivity or a combination of these. Another possibility is 
a change in the dust distribution, where a typical high- 
z LIRG/ULIRG may not be similar in structure to the 
strongly nuclear-concentrated local examples; instead the 
dust could be distributed over kpc scales, consistent with 
measurements of their extent in CO (e.g. Tacconi et al. 2006, 
Iono et al. 2009). 

Hints for evolution have also emerged through other 
studies of high redshift, IR-selected sources. For example, 
Rowan-Robinson et al. (2004, 2005) report results from the 
European Large Area ISO Survey (ELAIS) and the Spitzer 
Wide-area InfraRed Survey (SWIRE), which indicate that 
many IR-luminous galaxies in the 0.15-0.5 redshift range 
are best fitted with cool cirrus models. Examination of IR 
SEDs for z>0.5 ULIRGs, has revealed a higher cool to hot 
dust ratio than what is seen locally (Marcillac et al. 2006, 
Brand et al. 2008, Chakrabarti & McKee 2008). Sajina et 
al. (2006, 2008) present results based on mid and far-IR 
selected samples, which suggest that their sources are bet- 
ter described by cooler temperatures and high far/mid-IR 
ratios in contrast to local ULIRGs. Le Floch et al. (2005) 
find that ISOCAM-detected LIRGS at a median redshift 
of z~0.7 are characterised by half solar metallicities which 
could point to the IR SEDs being different from those of lo- 
cal LIRGs. This behaviour is an interesting complement to 
the evidence for evolution from mid-infrared spectroscopy 
and particularly the characteristics of Polycyclic Aromatic 
Hydrocarbon (PAH) features. The aromatic bands in high 
redshift galaxies are found to be similar in prominence and 
profile to the aromatic bands in local galaxies of lower lu- 
minosity (e.g. Pope et al. 2008, Rigby et al. 2008, Farrah et 
al. 2008), an important result as the PAH feature strengths 
in local galaxies are strongly correlated with emission from 
cold dust (Haas et al. 2002, Bendo et al. 2008). 

Predictions from several studies are also in agreement 
with the mounting evidence that the ISM, star-formation 
efficiency and dust properties of many high redshift in- 
frared galaxies significantly diverge from those in the local 
universe. Chemical evolution models (Calzetti & Heckman 
1999, Pei, Fall & Hauser 1999) predict a lower contribution 
in the UV SFR density and an increase in the average Ay 
by a factor of 2-3 at z~l-2, an outcome also supported by 
observations in Efstathiou & Rowan-Robinson (2003) and 
Le Floch et al. (2005). Le Floch et al. (2005) find a much 
stronger evolution in the infrared comoving energy density 
than in the UV over the same redshift range, evidence that 
dust reprocessing of optical light is more efficient at high 
redshift, but also tying in with an increase in star formation 
efficiency (SFR per unit cold gas mass) (Somerville, Primack 
and Faber 2001). 

The observations and results outlined above are likely 
to have considerable impact in accurately modelling the far- 
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Figure 7. Plot of the f35o/fs50 ratio versus redshift for LIRGs (top panel) and ULIRGs (lower panel). For reference the R08 tracks are 
also included (solid black curves). The flux density ratios corresponding to the 70fim sample are predicted using the SK07 models. Top 
panel: R08 LIRG tracks of logL /fl =ll, 11.25, 11.5, 11.75, all LIRG SMGs from Kovacs et al. (2006, K06) (large black diamonds) and 
Coppin et al. (2008, C08) (large black squares), all LIRGs in our sample at the original redshift z (blue filled circles), at a redshift z+0.7 
(blue open diamonds) and at a redshift z+1.5 (blue asterisks) — the change in redshift enables a comparison with the K06 and COS 
LIRGs. Lower panel: R08 ULIRG tracks of log L/fl=12, 12.25, 12.5, 12.75, all ULIRG SMGs from K06 (large black diamonds) and C08 
(large black squares), the ULIRGs in our sample at the original redshift z (red filled circles), at a redshift z+1.3 (red open diamonds) 
and at a redshift z+2 (red asterisks) — again the redshift change brings the 70 ^tm-detected ULIRGs to the redshift range of the K06 
and C08 ULIRGs. The horizontal error bars denote errors on photometric redshifts in the C08 sample. The vertical error bars for the 
70fj,m sample stem from the range of templates used to calculate f35o/f850 (see section f3.lt . Note that for a few SMGs, the uncertainty 
in Ljr could shift them to a higher luminosity class as some are very close to the upper luminosity limit. However, this does not alter 
the strong agreement in f35o/f850 colour between SMGs and the 70 fim sources, implying that we are detecting the lower redshift analogs 
of submm galaxies. 
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IR and submm background, as well as pinpointing the evolu- 
tion in the infrared luminosity function, where the cold SED 
component may play a more significant role than previously 
thought. We expect that cold galaxies such as those we iden- 
tify here, will be discovered in abundance by the Herschel 
Space Observatory (Pilbratt et al. 2004). The capabilities 
of Herschel's Spectral and Photometric Imaging Receiver 
(SPIRE, Griffin et al. 2007) and Photodetector Array Cam- 
era (PACS, Poglitsch et al. 2008) will provide detailed mea- 
surements in the 60-700 fira region. This will allow accurate 
mapping of previously unexplored territories such as the na- 
ture of the main contributors to the cosmic infrared back- 
ground and ISM properties at high redshift (Harwit 2004). 
The question as to whether a large fraction of IR-luminous 
galaxies at high redshift are cold, will explicitely be answered 
since the SED peak will be entirely defined at least up to red- 
shift z~3, enabling precise temperature measurements. In 
addition, with Herchel's capabilities the nature and proper- 
ties of dust at high redshift and its evolution through cosmic 
time will be accurately traced. 

5.2 The submm/IR connection 

SMGs, which are predominantly in the ULIRG and HyLIRG 
regime, are well-established as the most dust-rich galaxies in 
the Universe. Comparisons between the SEDs of SMGs and 
those of local and lower redshift (z<1.5) populations are 
therefore essential in the context of galaxy evolution, par- 
ticularly evident when modelling the behaviour of the IR 
and submm source counts, which necessarily requires the 
assumption of a SED shape. If SMGs are the high redshift 
analogs of local warm IR galaxies, the high submm counts 
on the faint end (with the limits of present surveys) could 
only be reproduced by evoking strong evolution well beyond 
z~3-4, as a warm SED moves the peak to shorter wave- 
lengths, reducing the submm flux for a given L/_r. This goes 
against the predictions of hierchical galaxy formation models 
which anticipate a downturn in space density of IR-luminous 
sources with Star Formation Rates (SFRs) >100M Q /yr by 
z~3-4. Instead, if submm emission arises from a significantly 
cooler SED, submm counts could be fit with galaxies of more 
moderate SFRs (Kaviani, Haehnelt and Kauffman 2003). 
Predictions by Efstathiou et al. (2000), Rowan- Robinson 
(2001) and Efstathiou &: Rowan-Robinson (2003) also sup- 
port this argument, suggesting that submm emission in high 
redshift galaxies could predominantly arise from cold dust 
and proposing that the quiescent cirrus component under- 
goes the same strong luminosity evolution as the starburst 
component. 

Up to now, studies of SMG SEDs converge on the idea 
that they are predominantly fit by cooler cirrus SED tem- 
plates (e.g. by Chapman et al. 2005, Kovacs et al. 2006, Pope 
et al. 2006, Coppin et al. 2008, Clements et al. 2008). Submm 
studies of high redshift ULIRGs and HyLIRGs, have shown 
that their intrinsic dust temperatures are distributed around 
a mean of T~35 K (Chapman et al. 2005, Kovacs et al. 2006, 
Pope et al. 2006, Coppin et al. 2008). Huynh et al. (2007) 
report very low typical dust temperatures (T=21-33 K) for 
a sample of L>5xlO n low-z SMGs. Interestingly, this cor- 
responds to the range in temperatures seen in the local SBs 
and LIRGs described in D00 and it is in contast to the higher 
temperatures observed in local ULIRGs — the Farrah et 



al. (2003) local ULIRG sample is fit by T~42K; see also 
Klaas et al. (1997) and Soifer et al. (1984) where the SEDs 
of Arp220 and NGC6240 are fit by modified blackbodies of 
temperature 47 K and 42 K respectively. Note that even if 
submm studies preferentially detect cold sources (e.g. Blain 
et al. 2004), many of these have total infrared luminosities 
comparable to the most energetic local galaxies, supporting 
the above evidence that a large number of high-redshift IR- 
luminous objects are colder than their local analogs. In fact 
the same conclusions are reached by Sajina et al. (2008) for 
a sample of high-redshift mid-IR selected sources, showing 
that this is not simply an artifact of the submm selection. 

Our results have established a connection between the 
far-IR and submm, showing that high redshift SMGs likely 
share similar properties with the lower redshift galaxies se- 
lected in Spitzer /MIPS surveys. We find that the predicted 
submm colours for the 70 /im-detected LIRGs and ULIRGs 
when redshifted to the 1.5<z<3.5 range, trace the same 
parameter space as the actual colours of submm-detected 
sources in the same luminosity and redshift range. This 
suggests that we have identified the missing link between 
SCUBA and Spitzer /MIPS galaxies, i.e. a population at in- 
termediate redshift (z ~ 1), detected by Spitzer but with the 
properties commonly seen in SCUBA galaxies. 



6 CONCLUSIONS 

We have examined the properties of a large sample of 
0.1<z<2 IR-luminous sources selected at 70 /im, correspond- 
ing to rest-frame 23-63 /im, a selection which should strongly 
favour starforming galaxies. Using available photometry in 
the 8-160/im range, we fit the Siebenmorgen & Krugel (2007, 
SK07) dust templates in order to calculate the total infrared 
luminosity and rest-frame wavelength of the SED peak, find- 
ing the former to be mainly in the 10 10 -10 14 Lq range and 
the latter at an average of ~105/im for the LIRGs and 
~92/im for the ULIRGs. As the main aim of this study 
is to compare our results with local sources of equivalent 
luminosity, we use data from the local IRAS Bright Galaxy 
Sample, again fitting the photometry with the SK07 mod- 
els. The locally-derived SED templates of Chary & Elbaz 
(2001, CE01) and Rieke et al. (2008, R08) arc also used as 
an additional comparative tool. We pay particular attention 
to any sources of bias that could affect our results and find 
nothing that could potentially predispose the outcome of our 
analysis. 

Our major findings are: 

• By examining the fi6o/f70 colour of the 70jj,m sample, 
we find that a large fraction of LIRGs and ULIRGs have 
higher fi6o/f70 than what is defined by the empirical CE01 
and R08 local galaxy SEDs, by at least a factor of 2 and up 
to 10 in extreme cases. As the fi6o/f70 continuum slope is 
directly related to the position of the SED peak wavelength, 
non-trivial differences between local and high redshift IR 
populations are implied. 

• For a given total infrared luminosity, a large fraction 
of the LIRG and ULIRG SEDs in the 70 ^m sample peak 
at wavelengths longer than 90 /jm, indicative of a signifi- 
cant cold-dust-associated far-IR excess not detected in local 
sources of comparable luminosity. This behaviour is directly 
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linked to evolution in dust and/or star- formation properties 
from the local to the high redshift universe. 

• We predict the f350 /fsso colours that the 70/^m-selected 
LIRGs and ULIRGs would have if found at the average red- 
shift of SMG populations (1 < z < 3) and find that they are 
on average equivalent to the observed SMG colours. This 
result suggests that we have potentially identified the z<l 
equivalents of the cold z>l sources discovered in blank field 
submm surveys. 

The strong evidence that we present here, drives the an- 
swer to the question that we defined at the start, i.e. whether 
there is a substantial number of infrared galaxies at high 
redshift with properties divergent from those of their local 
counterparts, to be yes. More specifically, our results have 
shown that these sources being colder rather than warmer. 
The emergent conclusion is that the evolution of the IR- 
luminous population with redshift is towards a larger dust 
mass or higher dust opacity and emissivity, a more extended 
dust distribution, higher star-formation efficiency or a com- 
bination of all. The imminent Herschel observational pro- 
grammes will be key in identifying and disentangling the 
main processes responsible. 
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